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We study by computer simulation a recently introduced generalised model of self-interacting 
self-avoiding trails on the square lattice that distinguishes two topologically different types of self- 
interaction: namely crossings where the trail passes across itself and collisions where the lattice path 
visits the same site without crossing. This model generalises the canonical interacting self-avoiding 
trail model of polymer collapse which has a strongly divergent specific heat at its transition point. 
We confirm the recent prediction that the asymmetry does not affect the universality class for a 
range of asymmetry. Certainly, where the weighting of collisions outweighs that of crossings this is 
well supported numerically. When crossings are weighted heavily relative to collisions the collapse 
transition reverts to the canonical ^-point-like behaviour found in interacting self-avoiding walks. 
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I. INTRODUCTION 

The collapse transition of a polymer in a dilute solution 
has been a continuing focus of study in lattice statistical 
mechanics for decades [TJ [5]. This transition describes 
the change in the scaling of the polymer with length that 
occurs as the temperature is lowered. At high tempera- 
tures the radius of gyration of polymer scales in a way 
swollen relative to a random walk: this is known as the 
excluded volume effect. At low temperatures a polymer 
condenses into dense, usually disordered, globule, with 
a much smaller radius of gyration. The interest in this 
phase transition has occurred both because of the motiva- 
tion of physical systems but also because of the study of 
integrable cases [31 [4] of lattice models have proved espe- 
cially fruitful in two dimensions. While the canonical lat- 
tice model of the configurations of a polymer in solution 
has been the self-avoiding walk (SAW), where a random 
walk on a lattice is not allowed to visit a lattice site more 
than once, an alternative has been to use bond-avoiding 
walks, or a self-avoiding trail. A self-avoiding trail (SAT) 
is a lattice walk configuration where the excluded vol- 
ume is obtained by preventing the walk from visiting the 
same bond, rather than the same site, more than once. 
The model of SAT was used initially to model polymers 
with loops 5 but has subsequently occurred in integrable 
loop models in two dimensions [I] . A model of collapsing 
polymers can be constructed starting from self-avoiding 
trails, known as interacting self-avoiding trails (IS AT). 
Here energies are associated with multiply-visited sites, 
and by favouring configurations with many such sites a 
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FIG. 1: An example of AISAT configuration with one crossing 
(m x — 1) and two collisions (m c = 2), associated with the 
Boltzmann weights t x and r c , respectively. The total number 
of doubly visited sites is m2 = 3. 



collapse transition can be initiated. 

Owczarek and Prellberg studied numerically the ISAT 
collapse on the square lattice by two different approaches 
[U [7] and in either case found a strong continuous tran- 
sition with specific heat exponent a — 0.81(3). Recently, 
on the triangular lattice Doukas et al. [5] found that by 
changing the weighting of doubly and triply visited sites a 
first-order transition can ensue or alternatively, depend- 
ing on the ratio of these weightings, a weaker second- 
order transition that mimics the collapse found in the 
canonical interacting self-avoiding walk (ISAW) model 
(also known as the 6*-point). They also found that the 
low temperature phase becomes fully dense rather than 
globular, depending on the choice of parameters. 

Recently, Foster [5] generalised the ISAT model on the 
square lattice by differentiating the type of doubly visited 
sites on the square lattice: a doubly visited site can be 
visited twice with the trail passing through itself, that is 
crossing, at the site, or alternatively as a result of two 
bends in the trail so that the trail "touches" or "collides" : 
see Fig. [T] 

We shall refer to this model as the asymmetric ISAT 
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model (AISAT). The study 9 using transfer matrices 
and the phenomenological renormalisation group of the 
AISAT predicted that the universality class of the sym- 
metric case extended to the asymmetric case. However, 
when crossings sufficiently dominate over collisions the 
results suggested the appearance of a first-order transi- 
tion. 

In this work we use Monte Carlo simulation to explore 
this AISAT model and the predictions of Foster [5] . We 
also explore the low temperature phase of the model and 
find that it is fully dense for a range of asymmetry in- 
cluding the symmetric case. 



II. ISAT 

The model of interacting trails on the square lattice 
is defined as follows. Consider the ensemble T n of self- 
avoiding trails (SAT) of length n, that is, of all lattice 
paths of n steps that can be formed on the square lattice 
such that they never visit the same bond more than once. 
Given a SAT ^ £ Tn, we associate an energy — e t with 
each doubly visited site and denote by mj(ip n ) the num- 
ber sites visited j times by ip n . We have n = mi + 2m,2- 
The probability of ip n is given by 



e /3e t m 2 (ip n ) 

z is -at ( T y 



(HI) 



where we define the Boltzmann weight uj t = exp((3s t ) 
and (3 is the inverse temperature l/fc^T. The partition 
function of the ISAT model is given by 



7 isat 



(T) 



The finite-length reduced free energy is 
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k„(T) = - log Z n (T) 



(II.2) 



(II.3) 



and the thermodynamic limit is obtained by taking the 
limit of large n, i.e., 



n(T) = lim K n (T). 



(II.4) 



It is expected that there is a collapse phase transition at 
a temperature T c characterized by a non-analyticity in 
«(T). 

The collapse transition can be characterized via a 
change in the scaling of the size of the polymer with tem- 
perature. It is expected that some measure of the size, 
such as the radius of gyration or the mean squared dis- 
tance of a monomer from the end points, i?^(T), scales 
at fixed temperature as 



R 2 n (T)~An 2 



(11.5) 



with some exponent v. At high temperatures the poly- 
mer is swollen and in two dimensions it is accepted that 



v = 3/4 [3J. At low temperatures the polymer becomes 
dense in space, though not space filling, and the exponent 
is v = 1/2. The collapse phase transition is expected to 
take place at some temperature T c . If the transition is 
second-order, the scaling at T c of the size is intermedi- 
ate between the high and low temperature forms. In the 
thermodynamic limit there is expected to be a singular- 
ity in the free energy, which can be seen in its second 
derivative (the specific heat). Denoting the (intensive) 
finite length specific heat per monomer by c„(T), the 
thermodynamic limit is given by the long length limit as 



C(T) = lim c„(T) . 



(II.6) 



One expects that the singular part of the specific heat 
behaves as 



C(T) ~ B\T C - T\ 



(II.7) 



where a < 1 for a second-order phase transition. The 
singular part of the thermodynamic limit internal energy 
behaves as 



U{T)~B\T C -T[ 



(II. 



if the transition is second-order, and there is a jump in 
the internal energy if the transition is first-order (an ef- 
fective value of a = 1). 

Moreover one expects crossover scaling forms [TU] to 
apply around this temperature, so that 

c n (T)~n a *C((T~T c )n*) (11.9) 

with < <j) < 1 if the transition is second-order, and 

c n (T)~nC((T~T c )n) (11.10) 

if the transition is first-order. From [10] we point out 
that the exponents a and <p are related via 

2-a=-. (11.11) 

Important for numerical estimation is the use of equa- 
tion (II. 9) at the peak value of the specific heat given by 



y 



peak 



T - T c )nV so that 



c P n ak (T) ~C peak n a4 > 



(11.12) 



where C peak = C(y peak ) is a constant. 

A previous study [7] of ISAT model on the square lat- 
tice has shown that there is a collapse transition with a 
strongly divergent specific heat, with 

a<p = 0.68(5) (11.13) 

and so the individual exponents have been estimated as 

4> = 0.84(3) and a = 0.81(3) . (11.14) 

At T = T c it was predicted [6] that 



R 2 n (T)~An (logn) 



(11.15) 
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III. ASYMMETRIC ISAT MODEL 

The Asymmetric ISAT (AISAT) model can be denned 
as follows. Consider the set of bond-avoiding paths T n as 
defined in the previous section. Given a SAT tp n G Tn, 
we associate an energy with each doubly visited site, as 
in ISAT, but we make a distinction between whether the 
trail crosses itself or not. We will call the former crossings 
and the latter collisions with associated energies — e x and 
— e c , respectively. For each configuration ip n £ T n we 
count the number m x (ip n ) of crossings and m c (ip n ) of 
collisions, see Fig. [T] Note that the total number of dou- 
bly visited sites is ni2 = m c +m x . We associate with each 
configuration a Boltzmann weight T ^ l!C ^' n ' T ^ lc ^' n ' where 
t x = exp(j3e x ), t c = cxp(/3e c ) and /3 is the inverse tem- 
perature 1/fceT. The partition function of the AISAT 
model is the given by 



Z n ( T x, T c) — ^ ' 



The probability of a configuration ip n is then 

p(ip n ;r x ,T c ) 



>x Tr. 



Z n { T x, T c) 

In line with Foster [pj let us define the variables 



and 



1 



(mi) 



(III.2) 



(III.3) 



(III.4) 



When we set t x = t c (x = 1, r = 1/2) the model reduces 
to the ISAT model, in which crossings and collisions are 
given the same weight. On the other end, if we set t x = 
(x = r = 0) configurations with crossings are excluded 
and our model reduces to the VISAW model |11) . 

The average of any quantity Q over the ensemble set 
of path T n is given generically by 



(Q){n;T x ,T c ) = ^ QWOpW 1 

i/>„6T„ 



n,T x ,T c ). (111.5) 



In particular, we can define the average number of cross- 
ings and collisions per site and their respective fluctua- 
tions as 



(III.6) 
(III.7) 



n 

K 

n 



n 

{m 2 c ) - (m c ) 2 



An important quantity for what follows is the proportion 
of the sites on the trail that are at lattice sites which are 
not doubly occupied: 



Pr> 



1 



((m c ) -I- (m x )) . 



(III.; 




FIG. 2: Density plot of the logarithm of the largest eigenvalue 
A mol of the matrix of second derivatives of the free energy 
with respect to t x and r c at length 500. 



Foster [9] predicted that the universality class of sym- 
metric ISAT at r = 1/2 extends to other values of r and 
further that there may be a change to a first-order tran- 
sition for large values of r. 



IV. NUMERICAL RESULTS 

We began by simulating the full two parameter space 
by using the flatPERM algorithm [12]. FlatPERM out- 
puts an estimate W n k of the total weight of the walks 
of length n at fixed values of some vector of quantities 
k = (fci, k2, ■ ■ ■ , ki). From the total weight one can access 
physical quantities over a broad range of temperatures 
through a simple weighted average, e.g. 



(O)n(T) 



(IV.l) 



The quantities kj may be any subset of the physical pa- 
rameters of the model. In our case we begin by using 
ki = irix and hi — m c . 

We have simulated AISAT using the full two-parameter 
flatPERM algorithm up to length n = 500, with 10 5 it- 
erations, collecting 1.2 • 10 10 samples at the maximum 
length. Following [TJ], we also measured the number 
of samples adjusted by the number of their independent 
growth steps ("effective samples") S e ^ ~ 3.9 • 10 s at 
the maximum length. To obtain a landscape of possi- 
ble phase transitions we plot the largest eigenvalue of 
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a(j) if the transition is second order. For the symmetric 
ISAT model, we estimated ad ss 0.64, which is a little 



T x = 1 



VISAW line 



Tc 



FIG. 3: Schematic diagram of the AISAT parameter space. 
The black dot on the diagonal line indicates the ISAT critical 
point. We have simulated along the dashed lines, as well as 
the VISAW line. 



the matrix of second derivatives of the free energy with 
respect to t x and r c at length n = 500 in Fig. [2] 

We notice that there is a strong peak in the fluctuations 
running in a line from r c 5 when t x — (r — 0) to 
t x fa 8 when t c = (r = 1). It is interesting to observe 
that the peak in the fluctuations becomes weaker as r 
increases from to I. 

Next, we have simulated four different one parameter 
slices of the AISAT model all up to length n = 1024. 
Their location, as indicated in Fig. |3J is as follows. 

• t x = (the VISAW model). With S ~ 1.8 -10 7 iter- 
ations, collecting 4.5 -10 10 samples at the maximum 
length (corresponding to S e ** ~ 5.2 • 10 8 ). 

• t x = 1 (the colliding model). With 7.8 • 10 6 itera- 
tions, collecting 2.6 • 10 10 samples at the maximum 
length (corresponding to S e ^ ~ 4.1 • 10 8 ). 

• t x = t c (the symmetric ISAT model). With S ~ 
4 • 10 6 iterations, collecting 7.6 • 10 9 samples at the 
maximum length (corresponding to S e *' ~ 5.8 • 
10 8 ). 

• t c = 1 (the crossing model). With S ~ 7.7- 10 6 iter- 
ations, collecting 2.7- 10 10 samples at the maximum 
length (corresponding to S eff ~ 2.7 • 10 8 ). 



A. Specfic heat 

We have begun by analysing the scaling of the spe- 
cific heat by calculating the location of its peak t% — 
argmax T c„(r) and thereby evaluating = c„(r^). In 
Fig. |4j we plot the peak values of the specific heat for the 
four models we have simulated. The exponent associated 



with the peak of the specific heat, see equation (11.12), is 



less that our previous estimate (equation (11.13)) based 
upon much longer length trails. We find for the VISAW 
(a<fi fa 0.69) and the colliding model (a(f> ~ 0.63) expo- 
nent estimates compatible with that of symmetric ISAT 
collapse at this length. This is a crucial result as it con- 
firms the prediction of Foster [9] that the universality 
class of ISAT extends from r = 1/2 to other values of r. 

We see that an attempted estimate of a<p for the cross- 
ing model over the full range of lengths is compromised 
by poor simulation results. Nevertheless there is clear 
curvature in the data and the estimate decreases with 
length. Using data from around length 400 onwards gives 
an estimate of 0.08. Given that this estimate would most 
likely decrease even further with increasing length, it is 
likely that the data is compatible with the behaviour of 
the #-point interacting SAW collapse transition; this has 
a negative value of the exponent a. 

Therefore our data indicates that the ISAT universality 
class extends to a range values of the parameter r mea- 
suring the asymmetry, as conjectured by Foster [5]. In- 
triguingly, rather than becoming first order for r near one 
as also conjectured in [S], the transition becomes weaker 
and potentially (9-like as r increases towards one. 



B. Study of the low temperature region 

To further investigate the proposition that changing 
the asymmetry over a range around r does not affect the 
nature of the collapse, except potentially near r near 1, 
we examine the low temperature phase. This approach 
proved fruitful in the case of triangular lattice extended 
ISAT [5] where the low temperature phase could be either 
globular or fully dense. The fully dense low temperature 
phase seems associated with either the ISAT universality 
class or a first-order transition, while the globular phase 
is associated with the much weaker (non-divergent spe- 
cific heat) ISAW (0-point) collapse transition universality 
class. Here we present evidence that for the square lat- 
tice symmetric ISAT, colliding and VISAW models the 
low-temperature phase is maximally dense, and that the 
density jumps discontinuously at the critical point. On 
the other hand for the crossing model the low tempera- 
ture phase seems not to be fully dense. This would be 
compatible with the conjecture that the collapse transi- 
tion for the crossing model is ^-point-like from a swollen 
phase at high temperatures to a globular phase at low 
temperatures. 

We have considered two different approaches to mea- 
suring the density. The first is indirect by measuring the 
proportion of sites of the lattice visited only once by the 
trail, and the second is using the radius of gyration to 
estimate the internal density of the polymer. Both lead 
to the same conclusions. 
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(a)The crossing model, r c = 1 



(b)The symmetric ISAT model, t c = t x 
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(c)The colliding model, 



(d)The VISAW model, 



FIG. 4: Double-logarithmic plots of the peak value of the specific heat against length for the four models defined at the 
beginning of Section 4. For three of the models the specific heat diverges strongly, while for the Crossing model the behaviour 
of the specific heat is markedly different: the data is consistent with a converging specific heat. 



1. Proportion of singly visited sites p n 

Following the analysis in [Sj we first measured the pro- 
portion p n of sites of the lattice visited only once by the 
trail. This provides a useful method for considering how 
dense our configurations are on average since an asymp- 
totic value of zero would imply that effectively all the 
sites occupied by the trail are doubly occupied. Dou- 
ble occupation of lattice sites implies that the surround- 
ing edges of the lattice are all occupied with bonds of 



the trail: hence the lattice is filled by the trail as the 
trail increases in length. At high temperatures it is easy 
to see that p n approaches a finite strictly positive value 
in the thermodynamic limit: this is connected with the 
swollen nature of the polymer as seen in the radius of 
gyration scaling. One would expect on physical grounds 
that Poo (T) would be a monotonically increasing function 
of temperature T. The question that arises is whether at 
low temperatures this value is zero and what kind of sin- 
gularity occurs in Poo{T) at the collapse transition. 
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FIG. 5: Plots of p n , the proportion of steps visiting the same site once, at a low temperature (on the left) and at the critical 
temperature (on the right). The scale n -1 ^ 2 chosen is the natural low temperature scale. The chosen temperatures for the plot 
on the left are r = 20 (crossings), r = 5 (IS AT), r = 7 (colliding), r = 8 (VISAW). At the critical temperature, the quantity 
p n is plotted against the crossover exponent obtained in [7] (1 — ct)<f> ~ 0.16. 



It is worth considering first what happens at the col- 
lapse point itself in the symmetric ISAT model. Thanks 
to the mapping between critical ISAT and Kinetic 
Growth Trails we know that critical limiting value for 
p n is exactly 1/5 0. The argument goes as follows: con- 
sider the case when a trail has formed a large n-step loop 
which occupies mi + to 2 = M lattice sites (this is al- 
ways the case as trails do not trap). Any site of this 
loop could have been the starting point. In order for this 
site to be visited only once, the loop must have closed 
at the first return visit, which occurs with a probability 
of 1/3. Therefore we find for large loops mi/M — > 1/3, 
from whence it follows p n = rtii/n —¥ 1/5. In Fig. [5] 
(a), we plotted p n against n -016 at the estimated criti- 
cal temperatures of our four sub-models. We have used 
(1 — ot)4> — 0.16 as the appropriate correction to scaling 
exponent from the estimated exponents for the symmet- 
ric ISAT model with a<f> w 0.68. For the symmetric ISAT, 
colliding and VISAW model this choice seems appropri- 
ate; in each case the estimated value of p^ is close to 
0.2. While the correction to scaling exponent may not 
be appropriate for the crossing model it is clear that any 
estimate of p^ is greater than 0.2. 

For low temperatures, as discussed above, if the trail 
fills the lattice asymptotically in a fully dense phase 
the portion of monomers of the trail not involved with 
doubly-visited sites of the lattice should tend to zero as 
n — » oo. In Fig. [H (b), we plotted p n against nT 1 ^ 2 
for our four sub-models. The plots suggest that for the 
symmetric ISAT, colliding and VISAW model, p n —> 



as n — > oo in the low-temperature region, which implies 
a maximally-dense phase: our extrapolation estimates 
for Poo have error bars encompassing zero. The cross- 
ing model seems to show a different behavior: our data 
suggest that p n tends to a non-zero value (around 0.25) 
These results form a clear picture in which the quan- 
tity Poo jumps discontinuously to zero at the critical tem- 
perature for the symmetric ISAT, colliding and VISAW 
models. In contrast, for the crossing model p^ remains 
non-zero for all temperatures. In Fig. [6j we have plotted 
the quantity p n at different length scales as a function 
of the temperature. The plots again suggest a common 
discontinuous behavior for the last three models, while 
the crossing model seems to show a continuous transition 
in line with the extrapolations described in the previous 
paragraph. 

These findings support the result of the specific heat 
analysis, namely that the universality classes of the sym- 
metric ISAT, colliding and VISAW models are likely to 
be the same, while the crossing model is clearly in a dis- 
tinct universality class. 



2. Density 

A more direct way of measuring of the density is to 
consider the quantity p — n/R 2 where R is the radius of 
gyration of the polymer. 

For any AISAT in the high temperature phase we ex- 
pect R 2 n ~ n 2v with v = 3/4, as for SAWs and SAT, and 
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(a)The crossing model, r c = 1 



(b)The symmetric ISAT model, t c = t x 
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FIG. 6: Density measured as the fraction p n of sites visited only once as function of the length scale n and of the temperature, 
for the four models defined at the beginning of Section 4. In (b) the dashed lines indicate the critical temperature and the 
corresponding value of p n . For three of the models p n tends to zero above a critical value of r, while for the Crossing model p n 
converges to a non-zero value for all values of r. 



therefore p n = n/R^ —> 0. From [S] we know that at crit- 
ical point of the symmetric model R n ~ n 1 / 2 log n, so 
the density p is zero also at the critical point. For all our 
models in the collapsed phase we expect that Poo{T) is 
non-zero at low temperatures. The natural question that 
arises is whether Poo(T) increases from zero at T = T c 
as the temperature is lowered in a continuous fashion, 
or whether it jumps discontinuously to a fixed maximum 
value at soon as the temperature is smaller than T c , being 



then constant for all T <T C . 

If we take our results for p n as a guide we would ex- 
pect that for the ISAT, the colliding and VISAW models 
Poo(T) the density would jump discontinuously on de- 
creasing T below T c while for the crossing model pco (T) 
would increases continuously from zero as T is lowered 
through T c . 

In Fig. [7] we plotted p n — njR\ as a function of the 
temperature. In the symmetric ISAT, colliding and VI- 
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(a)The crossing model, r c = 1 



(b)The symmetric ISAT model, t c = t x 
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FIG. 7: Plot of the density p n = n/R^ as function of the temperature, for the four models defined at the beginning of Section 
4. For three of the models there is a clear crossing point such that for small values of r the value of the density decreases with 
length, while for large values of r the value of the density increases with length. For the Crossing model no such crossing point 
exists, and the density decreases with length for all values of r. 



SAW models the density plotted for different lengths 
cross at a Boltzmann weight close to the expected col- 
lapse transition point, increasing to the right and de- 
creasing to the left. As argued above, the curves converge 
to zero for small Boltzmann weights (high temperatures) 
and are consistent with an increase to a temperature- 
independent constant for large Boltzmann weights (low 
temperatures). This is consistent with the behaviour of 
p n discussed previously. The crossing model seems to 



have a different behavior. Consistent with a varying lim- 
iting value of p n , the limiting value of the density p n 
does not approach a temperature-independent constant, 
either. 
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V. CONCLUSIONS 

We have studied a generalised model of interacting self- 
avoiding trails (ISAT) on the square lattice as proposed 
by Foster [5] where the weight associated with crossing- 
type interactions (t x ) and collision-type interactions (r c ) 
may differ. 

From our analysis, we can confirm the conjecture in 
[§] that the ISAT universality class extends over a region 
of asymmetry around t x — r c . We can conclude that 
this region extends down to t x = 0, which is also known 
as the VISAW model, and seems extends to some larger 
i~x > t c . Our simulation results for r c = 1 are compro- 
mised numerically by poor convergence: in fact at small 
lengths we see some evidence of multiply peaked proba- 
bility distributions but these seem to become unimodal 
at larger lengths. Importantly, the peak of the specific 
heat seems to decrease on increasing x = t x /t c . Indeed 
the specific heat of the t c = 1 model diverges with an 
effective exponent that is much smaller (albeit with large 
error bars) than the ISAT universality class would pre- 
dict and, in fact, the specific heat may not diverge: we 
would expect a more strongly diverging specific heat if 
a first-order transition occurs. Supporting these conclu- 
sions is our investigation of the low temperature phase 
for different asymmetry. For the symmetric ISAT and 
the VISAW model the low temperature phase seems to 



be fully dense as it is in the extended triangular lattice 
model in certain regimes. Also, in agreement with our 
tentative prediction for the t c — 1 case is the evidence 
that here the low temperature phase is no longer fully 
dense, which implies a globular ISAW-like low temper- 
ature phase. Putting this information together leads us 
to predict a ^-point-like collapse as occurs in interact- 
ing self-avoiding walks for x sufficiently large. Further 
numerical work is clearly needed to pin down the large 
x behaviour of the AISAT model, including where the 
change of universality class occurs. Finally, we point 
out that our evidence that the low temperature phase 
is fully dense for symmetric ISAT implies that there is a 
first-order characteristic of this transition, as predicted 
by Foster [5], even though there is no latent heat. 
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